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The coupling between superconductivity and othorhombic distortion is studied in vapor-grown 
FeSe single crystals using high-resolution thermal-expansion measurements. In contrast to the 
Bal22-based (Bal22) superconductors, we find that superconductivity does not reduce the or- 
thorhombicity below T c . Instead we find that superconductivity couples strongly to the in-plane 
area, which explains the large hydrostatic pressure effects. We discuss our results in light of the spin- 
nematic scenario and argue that FeSe has many features quite different from the typical Fe-based 
superconductors. 
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The interplay of structure, magnetism and supercon- 
ductivity has been a recurrent theme in the study of iron- 
based superconductors[T]. Among these systems, PbO- 
type /3-FeSe has the simplest crystallographic structure 
and a rich phase diagram|TJ. It undergoes a structural 
phase transition, similar to that of many parent com- 
pounds of the 1111 and 122 iron-based systems, at ~ 90 
K. In these later materials, structural and magnetic or- 
der track each other closely, which has led to the sugges- 
tion of a magnetic (spin-nematic) origin of the structural 
distortion^. In FeSe, however, no static magnetism is 
found at ambient pressure [3J H] and spin- fluctuations are 
found to be enhanced only at low temperatures^ , which 
raises the question of the origin of the structural transi- 
tion. FeSe becomes superconducting at a modest T c = 8 
KjB], yet the onset of superconductivity rises dramati- 
cally to ~ 37 K under hydrostatic pressure |S], which 
also induces static magnetic order|9 . Recently, super- 
conductivity even up to over 50 K was demonstrated in 
strained epitaxial thin films [10J. With its huge sensitiv- 
ity of T c to external pressure and the large separation 
between structurally distorted and magnetically ordered 
phase, FeSe is an intriguing system to study the phase 
interplay using pressure as a tuning parameter. Uniax- 
ial pressure effects, as can be studied by thermal expan- 
sion, are of special interest, because they are usually 
very anisotropic in the iron-based systems due to the 
layered crystal structures) 111 IT2"] . Further, these types 
of measurements provide a sensitive probe of the cou- 
pling between superconducting and orthorhombic order 
parameters (T3l [14] . 

Single-crystal growth of /3-FeSe, on the other hand, is 
complicated by the rich constitutional binary phase dia- 
gram of Fe and Se, where the tetragonal, superconducting 
/3-FeSe phase is located deep below the solidus iine[T5|. In 
consequence, crystal growth experiments from the melt 
or self flux result in hexagonal <5-FeSe, which undergoes 
a series of structural transformations and decomposition 
reactions on cooling to room temperature [TS]- Only if 



the initial Se content is low enough, can tetragonal j3- 
FeSe be obtained, however, only in the form of a platelet 
of hexagonal morphology which contains both, tetragonal 
/3-FeSe and magnetic Fe7Seg[15_. 

In this Letter, we report on high-resolution thermal- 
expansion of vapor-grown single-crystalline /3-FeSe. 
Growth directly in the tetragonal structure results in high 
quality single-crystals, which are well suited to study 
the anisotropy of pressure effects; a necessary comple- 
ment to hydrostatic-pressure studies in any non-cubic 
material. We show that the dramatic increase of T c 
under hydrostatic pressure arises from a reduction of 
the in-plane area, while T c is five times less sensitive 
to the c-axis length. Further, we demonstrate non- 
Fermi-liquid behavior in the low-temperature in-plane 
thermal-expansion coefficients, which slightly reduces the 
orthorhombic distortion and presumably arises from low- 
temperature spin-fluctuations. Surprisingly, orthorhom- 
bic distortion and superconductivity do not compete in 
FeSe, in contrast to underdoped 122 pnictides, which sug- 
gests that the structural transition may not be of mag- 
netic origin. 

Fe and Se powders were mixed in an atomic ratio 1.1:1 
and sealed in an evacuated Si02 ampoule together with 
a eutectic mixture af KC1 and AICI3. The ampoule was 
heated to 390° C on one end while the other end was kept 
at 240° C. After 28.5 days isometric FeSe crystals with 
tetragonal morphology were extracted at the colder end 
(Fig. [l](a,b)). At such low temperatures, samples form 
directly in the tetragonal state and do not undergo struc- 
tural transformations or decomposition reactions. Wave- 
length dispersive x-ray spectroscopy reveals an impurity 
level below 500 ppm, in particular there is no evidence 
for CI, Si, K or Al impurities. 

x-ray powder diffraction confirms the tetragonal struc- 
ture with lattice constants a — 3.7707(12) A and 
c = 5.521(3) A. Structural refinement with a 4-circle 
diffractometer using Mo-radiation yields a composition 
of Fe:Se=0.995(4):l (i. e. stoichiometric within the error 
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FIG. 1. (color online) (a,b) Photographs of tetragonal /3-FeSe 
grown using a low-temperature vapor-transport technique, 
(c) Temperature dependence of the magnetic susceptibility 
in a field of 1 T, applied parallel to the ab plane. The in- 
set shows the low-temperature data in a field of 20 mT. The 
screening is larger than —1 because of the demagnetization 
effect. 



bar) and a structural z parameter of z = 0.26668(9). No 
indications for interstitial atoms were found. Fig. [I] (c) 
shows the temperature dependence of the magnetic sus- 
ceptibility with field applied parallel to the ab plane mea- 
sured in a vibrating sample magnetometer. A small but 
sharp kink, which we associate with the structural tran- 
sition, is observed at 87 K. The superconducting transi- 
tion has a sharp onset at 8 K and is broadened by the 
relatively high applied field. High-resolution thermal ex- 
pansion was measured in a capacitance dilatometer|16|. 
in which the sample is pressed against one plate of a 
plate-type capacitor with a force of ~ 0.2 N, directed 
along the measured sample length. As for underdoped 
BaFe9As9 |13l[T3] . this force can be used to in-situ detwin 
samples of FeSe below their tetragonal-to-orthorhombic 
phase transition (from space group P4/nmm to Cmma). 
In particular, if the thermal expansion along [110]t 
(tetragonal notation) is measured, samples are detwinned 
by the applied force and the (shorter) orthorhombic a 
axis is measured[17 . Measuring thermal expansion along 
the tetragonal [100]t direction yields, ideally, the aver- 
age of a and b axis. Under this assumption, the thermal 
expansion of the orthorhombic b axis can be inferred [14"]. 

Fig. [2] shows the thus obtained relative sample-length 
changes — L ^ T L~(3oiyK) ^ ano - uniaxial thermal- 

expansion coefficients oij = jr^f, where the index i 
stands for the direction. The data are in good 

agreement with previous x-ray studies[18 and show clear 
evidence for a second-order tetragonal-to-orthorhombic 
phase transition at T s = 87 K. We note that the 
c-axis anomaly at T s is unusually small when com- 
pared to underdoped Ba(Fe,Co) 2 As 2 (Co-Bal22)|T3] or 
BaFe 2 (As,P) 2 (P-Bal22)[Tl]. The very small anomaly in 
the volume average of the thermal expansion indicates 
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FIG. 2. (color online) (a) Relative length changes along the 
three measured directions (orthorhombic a axis, in-plane aver- 
age and c axis) and inferred fo-axis length change (continuous 
lines). For comparison, the corresponding data from x-ray 
diffraction|18] are given (circles), (b), (c) and (d) show the 
6-axis, in-plane average and a-axis length change close to T c , 
respectively, on a magnified scale, (e) Uniaxial thermal ex- 
pansion coefficients of FeSe along three measured directions, 
with the thermal-expansion along the 6 axis and the volume 
average inferred from the measurements. The inset shows the 
data close to T c on a magnified scale. 



that T s does not couple strongly to hydrostatic pressure. 
No distinct second, potentially magnetic, phase transi- 
tion is observed below T s . At T c = 7.75 K, the dis- 
continuity in the in-plane thermal-expansion coefficients 
Ac*i (kink in ALj(T)) clearly confirms a sharp, bulk su- 
perconducting transition. Acvi is related to the uniaxial 
pressure derivative of T c via the Ehrenfest relationship 



dT c _ VmAOj 

dpi ACp/T c 



(1) 



V m = 23.34 cm 3 /mol is the molar volume and AC P /T C — 
9.45(30) mJ mol _1 K -2 is the specific heat jump, which 
we take from Reference [HJJ We thus obtain d ^ a = 

2.6(3)K/GPa for the in-plane average, g = 2.2(5) 
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K/GPa, dT < 



- 3.1(1.1) K/GPa and ^ = 0(0.5)K/GPa. 
The hydrostatic pressure derivative of T c is simply given 
by the sum of the uniaxial components = |^ + 

2^ + = 5.3(1.2) K GPa -1 and is in good agreement 
with the initial slope of direct measurements, which yield 
g = 6-7K GPa _1 [9l [20]. It is clear from our results 

that the comparatively large ^ in FeSe arises not from 
particularly large uniaxial components but from lack of 
their cancellation. The in-plane derivatives are compara- 
ble in size to slightly overdoped Co-Bal22 |lll I12| . How- 
ever, in Co-Bal22, in-plane and c-axis pressure deriva- 
tives have opposite signs and largely cancel in the hydro- 
static average (with the negative p c derivative slightly 
prevailing) |11[ I13j . while the p c derivative is approxi- 
mately zero in FeSe. The structural tuning parameter of 
the 122-systems is the c/a ratio|i T2TfT3] . In FeSe instead, 
the in-plane distance alone appears to be the tuning pa- 
rameter, which couples strongly to hydrostatic pressure. 

Basically the same picture emerges when considering 
the uniaxial strain derivatives of T c , = Si^i^J) 
which can be calculated if the set of elastic constants cy 
= 1 — 3) is known. For tetragonal FeSe, these con- 
stants have been calculated using DFT|5T] to en = 95.2 
GPa, ci2 = 48.8 GPa, c 13 = 13.9 GPa and c 33 = 39.5 
GPa|22 . Using the above values (allowing for an er- 



ror of 10%) we find 



-399(177) K. 



-358(142) 



K and ^ = -73(49) K. It is evident that T c depends 
sensitively on the in-plane lengths. The corresponding 
uniaxial Griineisen parameter is Tc = —49(12), i. e. 
shrinking (a + b)/2 by 1% increases T c by ~ 50%. Note 
that a and b decrease by ~ 1.5% between ambient pres- 
sure and 7 GPa, where the highest T c is reached|8 . T c 
is five times less sensitive to changes of the c-axis length 
( d ^" Tc = -9(6)). It is striking that the c-axis length has 
such a small effect on T c , especially since the Se height 
(i.e. the c-axis length times the internal z-parameter) was 
found to correlate closely with T c [23]. The z-parameter 
however, may depend in a complicated manner both on 
in-plane and c-axis lengths. For a detailed investigation 
of its relation to uniaxial pressure effects, this dependence 
would have to be established. 

Fig. |3jpresents our results concerning the interplay be- 
tween orthorhombicity, magnetic fluctuations and super- 
conductivity in FeSe. Fig. [3] (a) shows the orthorhombic 
order parameter S = \a — b\/(a + b) s» \a — b\/2a,Q of FeSe, 
undoped Bal22 and underdoped Co-Bal22 (4.5% Co- 
content), all computed from our thermal-expansion mea- 
surements. Although no magnetic transition is found, 
6(T = 0) of FeSe is of similar magnitude as in the 122- 
systems. While S(T) is reduced below T c in the 4.5%Co- 
Bal22, there is no discernable feature in S(T) of FeSe at 
T c (Fig. [3](b,c)). Curiously, however, 5{T) has a weak 
maximum at ~ 12 K, which we will discuss before we 
address the response of S(T) at T c in more detail. 

The maximum of S(T) is caused by the sign change of 
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FIG. 3. (color online) (a) Orthorhombic order parameter 5 = 
fe— a\/2ao of FeSe (orange line), undoped Bal22 (dashed black 
line) and underdoped Co-Bal22 (black line) computed from 
thermal-expansion measurements, (b) and (c) show 5(T) of 
FeSe and Co-Bal22, respectively, below 20 K on a magnified 
scale. Note that the scale in (b) is ten times smaller than 
in (c). (d) a a ,t/T of FeSe (red and purple lines) and a a /T 
of underdoped P-Bal22 (dark green line). For clarity, data 
below T c are shown in a lighter color. The arrow points at 
the crossing of the normal-state in-plane thermal-expansion 
coefficients of FeSe. (e) and (f) show ^ of FeSe and Co- 
Bal22, respectively, close to T c . 



both a a (T) and a b (T) at - 12 K (see inset Fig. ge)), 
which points to an additional low-temperature contribu- 
tion to the thermal expansion. This contribution, which 
appears to be diverging down to T c , becomes evident in 
a plot of <y ^- (Fig. [3] (d)). For a Fermi-liquid, one ex- 
pects a constant ^-term at low temperatures, which is 
directly related to the uniaxial pressure derivative of the 
Sommerfeld coefficient, as seen for P-Bal22 with 20% 
P-content[14] (Fig. ^d)). The non-Fermi-liquid charac- 
ter of the thermal-expansion of FeSe becomes apparent 
in this comparison. A new energy scale with a negative 
(positive) contribution to ^ (^) emerges below ~ 30 K 
and causes these coefficients to diverge. Note that nei- 
ther the Knight-shift, nor specific heat have previously 
shown indications of non- Fermi- liquid behavior [5J IT!?] . 
However, thermal-expansion is expected to be an espe- 
cially sensitive probe for locating such non-Fermi-liquid 
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behavior[2U and enhanced low-temperature spin fluctu- 
ations have reavealed that FeSe is on the brink of a mag- 
netic phase transition near T = (5]. Our observation 
of non-Fermi-liquid behavior is probably related to this 
quantum critical point and the new energy scale may 
be linked to the low-temperature spin fluctuations. It 
is then, however, curious that these fluctuations cause a 
reduction of the orthorhombic order parameter of FeSe. 
In the 122-pnictides, in contrast, the onset of magnetism 
enhances <$(T)[25], suggesting that the relation between 
magnetism and orthorhombicity may be different in FeSe. 

The very small et6-plane anisotropy of Aa^ of FeSe at 
T c implies that 5(T) has no more than a tiny anomaly at 
T c , as can indeed be seen in Fig. [3](b). Only a small kink 
in the temperature derivative (Fig. J3^d)) is observed, 
which is such that i5 will have a tendency to be somewhat 
larger in the superconducting state than in the normal 
state. In particular, there is no indication of competition 
between orthorhombicity and superconductivity in FeSe. 
Such a competitive coupling was clearly observed in un- 
derdoped Co-Bal22 (see Fig. [3^c) and Ref. l2"o| and other 
doped Bal22-compounds|27 and results in a reduction of 
8 below T c , or equivalently in a positive anomaly in 4= 
(see Fig. [3jf ) ) . This coupling was proposed to arise, ul- 
timately, from the competition between magnetism and 
superconductivity [25]. In this scenario, magnetic spin 
fluctuations, which give rise to "nematic" order that, in 
turn, induces an orthorhombic distortion via magneto- 
elastic coupling, are weakened by the onset of supercon- 
ductivity. It is an important question which part of this 
scenario is not valid for FeSe. Our results suggest that 
either superconductivity does not interact strongly with 
spin fluctuations or that the structural transition has a 
non-magnetic origin. 

In conclusion, high-resolution thermal-expansion mea- 
surements of vapor-grown /3-FeSe have revealed a number 
of unusual properties, which suggest that FeSe is not a 
typical iron-based superconductor. The structural tun- 
ing parameter of FeSe is the in-plane area, compared to 
the c/a ratio in 122-systems, which explains why T c cou- 
ples so strongly to hydrostatic pressure. In underdoped 
122 pnictides, magnetism and orthorhombicity are coop- 
erative and compete with superconductivity. In FeSe, on 
the other hand, T c and magnetism both increase under 
hydrostatic pressure, suggesting that they act coopera- 
tively. Further, superconductivity does not compete with 
orthorhombicity, which raises the question of the origin of 
the orthorhombic phase transition. Also, we show that a 
new energy scale, presumably associated with spin fluctu- 
ations, emerges at low temperature and slightly reduces 
the orthorhombic distortion below 12 K. Further stud- 
ies on the nature of these different phases, which exhibit 
such an unusual interplay, will be of great interest. 
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